The Universal Negative Poisson Ratio of Self- Avoiding Fixed-Connectivity 

Membranes 
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We determine the Poisson ratio of self-avoiding fixed- 
connectivity membranes, modeled as impenetrable plaque- 
ttes, to be a = —0.37(6), in statistical agreement with 
the Poisson ratio of phantom fixed-connectivity membranes 
a = —0.32(4). Together with the equality of critical expo- 
nents, this result implies a unique universality class for fixed- 
connectivity membranes. Our findings thus establish that 
physical fixed-connectivity membranes provide a wide class 
of auxetic (negative Poisson ratio) materials with significant 
potential applications in materials science. 



Fixed-connectivity (also known as polymerized, teth- 
ered or crystalline) membranes are fluctuating and flex- 
ible fishnet-like two-dimensional surfaces with nodes of 
fixed coordination number (for two recent reviews see 
PIP))- Physical examples include such naturally occur- 
ring structures as polymerized Langmuir-Blodget films 
[pLE] and the spectrin/ actin cytoskeleton Q of erythro- 
cytes (mammalian red blood cells). A wide variety 
of additional examples is discussed in 0. Current 
advances in soft condensed matter experimental tech- 
niques suggest the likelihood of many new realizations of 
fixed-connectivity membranes such as cross-linked DNA 
networks as well as composite structures that include 
fixed-connectivity membranes as fundamental ingredi- 
ents. One universal and remarkable feature of the low- 
temperature (so-called flat) phase of non-self-avoiding 
(phantom) fixed-connectivity membranes is that they ex- 
pand transversely when stretched longitudinally . In 
other words, they exhibit a negative Poisson ratio Jl0| . 
Such materials have been dubbed auxetic jO]. 

In this letter we estimate, via Monte Carlo simula- 
tions, the Poisson ratio of physical self-avoiding fixed- 
connectivity membranes. We establish that they are also 
auxetic materials with a Poisson ratio and roughness ex- 
ponent in statistical agreement with those of flat phan- 
tom membranes. Thus there appears to be a unique 
universality class of flat fixed-connectivity membranes, 
whether they arise from high bending rigidity or self- 
avoidance. Direct experimental measurements should be 
able to measure this negative Poisson ratio and test uni- 
versality. The remarkable properties of auxetic mem- 
branes suggest a rich new avenue of exploration in mate- 



rials synthesis. 

A fixed-connectivity membrane may be modeled as an 
elastic surface with bending rigidity and self-avoidance 
|l2| [l5[ . The free energy is composed of an elastic and a 
self-avoiding contribution: 



Ffcm F e i -\- F sa . 

The elastic free-energy F e i is given by 
F el {u,h) = I d 2 x 



k A 
-(d a d p h) 2 + fiu 2 a/3 + -u 2 ai 



(1) 



(2) 



where u denotes phonon modes, h denotes height modes, 
k is the bending rigidity and A and fi are the classical 
Lame coefficients. The self-avoiding free-energy we take 
to be of the Edwards form 



F sa = ~ I d 2 x 



d 2 y<5 3 (f(x)-r(y)) , 



(3) 



where b determines the strength of self-avoidance. The 
strain tensor u a p fT(| is related to the embedding f*(x), 
defining the membrane by 



r(x) = x + u(x) + zh(x) 
u a p = -^{daUp + d/3U a + d a hdph) 



(4) 



Combined efforts both on the analytical and numerical 
side have led to a complete clarification of the phase dia- 
gram of the phantom case (b — 0) and detailed estimates 
for the critical exponents, as shown in Fig. [l] and Table |. 
The phase diagram consists of a crumpled phase (associ- 
ated with the Gaussian fixed point GF) and a flat phase 
(associated with a flat phase fixed-point FL), with an in- 
termediate infrared- unstable crumpling transition (CT), 
as depicted in Fig. |]. 

The self-avoiding model of Eq. ([!]) with no bending 
rigidity (k = 0) has proven to be tractable numerically 
]l6| . The model possesses a unique infrared fixed point 
(SA) describing a flat phase, with detailed results for 
critical exponents given in Table [j]. From the analytical 
standpoint there is evidence for a unique S A jl7],[l8| , but 
the calculations are inconclusive on whether it describes 
a flat or a crumpled phase fll9|-p4j . 

A look at Table || reveals that the critical exponents 
and the Poisson ratio of both phantom and self-avoiding 
membranes coincide within the error bars quoted, im- 
plying that the FL and the SA are equivalent. This is 
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TABLE I. Critical exponents and Poisson ratio of fiat 
fixed-connectivity membranes in the phantom (MC:Monte 
Carlo; SCSA:Self-consistent screening approximation) and 
the self-avoiding case (MC-BS:Monte Carlo with Balls and 
Springs models; MC-IP:Monte Carlo with impenetrable mod- 
els). 



FIG. 1. The renormalization group flows in the 
two-dimensional space of couplings for a fixed-connectivity 
membrane with bending rigidity (ft) and self-avoidance (6). 
The phantom model (b — 0) has two infrared stable fixed 
points, the crumpled phase (GF) and the flat phase (FL), 
with an intermediate continuous crumpling transition as- 
sociated with the infrared-unstable fixed point (CT). The 
pure self-avoiding model with no microscopic bending rigidity 
(/t = 0) has a infrared stable self-avoiding fixed point (SA). 
There is a line of equivalent fixed-points joining the FL and 
the SA (the solid (red) line), thus defining a redundant direc- 
tion in ft — b space. 



a very surprising result since it means that the same 
long wavelength limit is reached via two very different 
routes: either sufficiently large bending rigidity or strong 
self- avoidance in the absence of bending rigidity. This 
is even more remarkable considering the very different 
short-distance structure of the two models: a membrane 
with large bending rigidity is very smooth at short dis- 
tances while a purely self-avoiding membrane in the ab- 
sence of bending rigidity is extremely rough. 

In light of our results, we suggest the renormalization 
group (RG) flows depicted in Fig. [l]. Given these flows 
and the equivalence of the SA and the FL, there must be a 
full line of equivalent fixed points 6(ft) joining them, with 
no RG flow on the line (corresponding to a marginal di- 
rection). Since these are the infrared stable fixed points 
of the system, all the relevant physics is described by 
this line. The crumpling transition present for phantom 
membranes may be reached only by an extremely pre- 
cise tuning of the parameters involved in the problem. 
This would make the crumpling transition in the model 
described very difficult to verify experimentally. 

We see that the combination of fixed-connectivity (in- 
tegrity of the lattice) and self-avoidance sufficiently re- 



stricts the entropy of crumpled configurations as to de- 
stroy the crumpled phase. It was already observed in |2S| ] 
that next-to-nearest neighbor self-avoidance, discretized 
by hard-sphere potentials, induces a positive bending 
rigidity. On the other hand the impenetrable plaquette 
model treated here is very flexible, since only strictly self- 
intersecting configurations are forbidden, and hence the 
physical mechanism flattening the membrane is clearly 
more general than the simple induced bending rigidity 
discussed above. 

We performed a Monte Carlo simulation of a suitable 
discrete version of the model (Eq. (|l|)), as described in 
|l6| . The Poisson ratio a for a two-dimensional system 
deformed from its mean length / by SI is determined by 



Sw/w 
51/1 



(5) 



where Sw/w measures the fractional change in the trans- 
verse extent (width) of the system. In || it is shown that 
linear response theory gives 



(Uyy }c 



(6) 



where (u 2 ) c is the connected statistical average over 
Monte Carlo configurations and 71 is the spatial average 
over the surface. 

Our results are presented in Fig. ||(A) and show that 
the Poisson ratio a = —1/3 obtained in simulations of 
flat phantom fixed-connectivity membranes is also a good 
estimate of the Poisson ratio for the pure self-avoiding 
model. Error bars have been calculated with the Jacknife 
method 

To calibrate the influence of the boundary on our re- 
sults we calculated the Poisson ratio excluding concentric 
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FIG. 2. Poisson ratio of a self-avoiding fixed-connectivity 
membrane as a function of system size (A). Poisson ratio 
of reduced lattices compared with non-truncated ones (B). 
The straight dashed line indicates in both cases the SCSA 
analytical result a = —1/3. 



outer shells of nodes neighboring the boundary. In par- 
ticular we performed the numerical analysis discarding 
shells of boundary nodes of increasing extent. The re- 
sults are presented in Fig. ||(B). 

The exclusion of larger shells of boundary sites slightly 
increases the absolute value of the Poisson ratio. For suf- 
ficiently small reduced linear size, |er| begins to decrease 
due to finite size effects. There is consequently competi- 
tion between boundary and finite size effects. Finite size 
effects become important only when a sizeable fraction 
of nodes near the boundary are excluded. Note that we 
can also make a consistent check of our analysis by sys- 
tematically excluding shells of nodes from the boundary 
in towards the center of the lattice and comparing the 
lattices of reduced size thus obtained with equal volume 
non-truncated lattices. 

Thus we can compare, for example, the L = 65 result 
with that from the reduced L = 95 lattice and likewise 
the L = 49 result with that from the reduced L = 95 and 
L = 65 lattices. The matching given by this comparison 
is consistent. In fact we are able to reproduce the L = 
65, L = 49 and L = 33 results simply by reducing the 
L — 95 lattice. The deviations found in this comparison 
are another measure of boundary effects. 

Traditional materials get thinner when stretched and 
fatter when squashed, since it is typically difficult to in- 





FIG. 3. Mechanical model of an auxetic material: (a) in 
the absence of applied stress and (b) under applied lateral 
stress T. The lateral stretching accompanying the applied 
stress forces the material out in the transverse dimension. 



crease their volume very much when deformed In the 
unusual world of auxetics the opposite happens, with a 
number of interesting implications and potential appli- 
cations. There are several well known auxetic materials. 
The earliest example, dating from more than a century 
ago, is that of a pyrite crystal J32], which has a Poisson 
ratio in certain directions of a ~ —0.14. More recently 
some polyester foams under certain pressure conditions 
have proved to be isotropic auxetic materials with Pois- 
son ratios as large as a ~ —0.7 p3| . A nice mechanical 
model for auxetic materials was given in [|ll| (see Fig. |^) . 
One of the rare naturally occurring auxetics is Si02 in 
its a-crystobalite phase P^ , P5| . 

The underlying mechanism driving fixed-connectivity 
membranes auxetic has some similarities to that illus- 
trated in Fig. |§|. Submitting a membrane to tension will 
suppress its out-of-plane fluctuations, forcing it entropi- 
cally to expand in both in-plane directions. More physi- 
cally, the out-of-plane undulations renormalize the elastic 
constants (the Lame coefficients), in such a way that the 
long- wavelength bulk modulus is less than the shear mod- 
ulus, which is the signature of a two-dimensional auxetic 
material. 

Auxetic materials have desirable mechanical properties 
such as higher in-plane indentation resistance, transverse 
shear modulus and bending stiffness. They have clear 
applications as sealants, gaskets and fasteners. They may 
also be promising materials for artificial arteries, since 
they can expand to accomodate sudden increases in blood 



*A volume preserving deformation has Poisson ratio 0.5. 
Any value less than 0.5 involves some increase in volume un- 
der deformation. A negative Poisson ratio implies a very large 
volume increase. 
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flow. 

We can model a realistic fixed-connectivity membrane 
with an elastic free energy and either large bending rigid- 
ity or self-avoidance. This is of practical importance in 
modeling since, for example, we may replace the more 
complicated non-local self-avoidance term with a large 
bending rigidity. It remains an important theoretical 
challenge to verify this conclusion analytically. 

In this letter we have ignored the role of topological 
defects. We think that a sufficiently large defect density 
may affect the actual value of the Poisson ratio, but a 
detailed discussion of this topic is beyond the scope of 
this paper. 

It is our hope that the results presented in this letter 
will encourage materials scientists and condensed matter 
experimentalists to further study the elastic and mechan- 
ical properties of fixed-connectivity membranes and, in 
particular, to measure the Poisson ratio of these novel 
systems. 
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